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PREFACE 


ASTER Consulting Assoointos, Ino, A?AS Marietta Drive, BlnRhamton, 

New York, was awarded a eoutraet (NAS 9-16505), Involving amu-oxlinately 
1,5 months of effort, to review the crop canopy models. 

The technical requirements of the cimtract were: 

(1) to objectively review the oxlstinp, crop canopy models and their 
extensions Cor physical, biophysical and mathematical assiuuptlons, 
visor, and validity. 

(2) to review any other models that are of relevance to this area, 

(3) to include sensitivity of results to assumptions and extendvahlUty 
of these models for realistic conditions. 

(A) to suRRVst areas for hnprovement and the need for new models. 

The author was to bo tlv. primary reviewer , This is the final report 
on the contract . 
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ABSTRACT 


Various models for calculating crop canopy reflectance, in the 
visible and infrared wavelengths, from the optical and geometrical 
properties >f a canopy and its constituents have been reviewed. 

The review includes a discussion of radiative transfer equation, and 
both analytical and numerical crop reflectance models which are 
manifestations of the solution of this equation. Recommendations are 
made for further work in modeling of canopy reflectance. These 
recommendations include: 

. extensive testing of various models using a common data base 

. development of a simple and still fairly realistic crop reflectance 
model, involving adaption of existing models, as a short term effort 

. development of a more realistic and numerical model as a long term 
effort 

. development of a canopy reflectance model with time implicit i:i it by 

combination of a canopy reflectance model with a crop growth model. 

. Investigation of the potential of various canopy reflectance models 

in determining canopy parameters of importance from reflectance data, 
i.e. model inversion 

. adaption of other invertible reflectance models thin layer 

system model to crop identification problem. 

• detailed study of other modeling areas relevant to crop canopy reflec- 
tance, including modeling of reflectance from single canopy component, 

modeling of reflectance in thermal and microwave regions and modeling 

atmosphere, to determine the strategy for future work. 
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I. INTRODUCTION 


Wlien daylight is Incident on a vegetation (crop) canopy, it is 
scattered and reflected, and its direction and spectral composition 
are altered in a complex manner by the vegetation.' Part of this 
altered radiation is remotely sensed by Landstat. It is hoped that 
with multispectral measurements one can determine crop canopy para- 
meters which can assist in 
. crop identification 
. crop growth stage determination 
. crop quality or abnormality determination 
and 

. eventually crop yield calculations. 

The role of crop canopy reflectance modeling towards achieving this 
objective is stated rather well by Suits (1981), wlio proposed a crop 
canopy model which has been developed and investigated by many scientists 
and is widely used. According to him; 

' A canopy reflectance model provides the logical connection between 
the botanical features of the canopy, the geometry of the radiometric 
interaction and the resulting alteration in the reflected radiation. Such 
a model allows one to understand the reasons for the alterations and to 
calculate the magnitude and trends of these alterations caused by the 
botanical features and the geometry of the interaction. The validity of 
inferences as to the existence of important agronomic features from the 
detected altered radiation may be tested on theoretical grounds'. 
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Ovsr the last 15 years or so, several canopy reflectance models 
have been developed. These models represent either an approximate or 
a numerical attempt to solve what is known as radiative transfer equation 
whicl is a macroscopic manifestation of the interaction of radiant 
energy with matter. In Section II, we provide a general perspective of 
the interaction of electromagnetic waves with crop canopies including 
radiative transfer theory. In Section III, we summarize those crop 
reflectance models which are based on approximate but analytical solution 
of the radiative transfer equation. We shall refer to them as analytical 
models. In Section IV, we summarize numerical crop reflectance models where 
an attempt is made to numerically solve the radiative transfer equation. In 
the discussion of both types of models , vje v’ill emphasize the key assumptions 
made, and point out agronomic variables which are used in the calcuation of 
reflectance, with the hope that it will assist the reader in assessing the 
usefulness of a model in estimating agronomic variables from the reflectance 
data. For details of various models, the readers are referred to the original 
papers. The readers may also find an unpublished report ’MRS Literature Survey 
of Bidirectional Reflectance and Atmospheric Corrections II. Bidirectional 
Reflectance Studies Literature Review' prepared for NASA-Goddard Space Flight 
Center, Greenbelt, MD. by J.A. Smith and K.J. Ranson useful. This extensive 
report (about 200 pages) dated August 1979 provides an excellent comprehensive 
review of previous work in scene bidirectional reflectance, an extensive 
bibliography and abstracts of key papers. 

I 

In Section V, we will recommend a strategy for further work in the areas 

of crop reflectance modeling. This strategy includes comparative testing of 

various models, initiation of work involving inverting the models, i.e. to 
develop procedures for obtaining agronomic variables from the reflectance 

data using these models, and the desirability of new models. 


II. 


INTERACTION OF ELECTROMAGNETIC WAVES WITH CROP CANOPIES - GENERAL PERSPECTIVES 


The theoretical basis for understanding the interaction between electro- 
magnetic radiations and crop canopies is the radiative transfer theory, also 
called transport theory. It is a macroscopic theory of the interaction of 
radiant energy with matter. It describes the observed phenomena of light 
scattering, absorption, and polarization effects, but without regard to the 
classical electromagnetic theory. This theory has been developed, and applied 
extensively by astrophysicists, earth and atmospheric scientists for studying 
steller or planetary atmosphere, earth surfaces and oceans. One of the classic 
and encyclopedic work in the field is due to a well known astrophysicist, 

S. Chandrasekhar (1950) . Since then many texts and monographs have appeared 
emphasizing different aspects of the theory, including a more recent and 
rather comprehensive set of two volumes by Ishimaru (1978 a, b). 

The mathematical apparatus of the radiative transfer theory, though 
conceptually straightforward, can not be easily applied to the vegetation case 
because cf a number of unusual features of the vegetation. To appreciate these, 
from the pedagogical point of view, it is desirable to provide a quick review 
of the radiative transfer theory, adapted from Ishimaru (1978a) and to a lesser 
extent from Smith and Ranson (1979) . 

The transport theory describes the propagation of intensities in randomly 

distributed particles in terms of specific intensity I(r,s), which is in general 

function of position r and direction s in a three-dimensional space. Specific 

intensity, also called radiance or brightness , is the average power flux density 

within a unit frequency band centered at frequency v within a unit solid angle 

—2 “1 “1 2 

and is measured in W m~ Sr H (Watts/meter /steradian or solid angle/Hertz) . 
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This quantity satisfies the so called equation of transfer 


. _ f 3- T Cr. 's )+ p.r 3-Cr/O 6 (. r, S J 

A 5 ^ (2.1) 

In this equation, p is the number of particles per unit volume with which 
the incident radiation interacts and a is the total of scattering and absorption 
cross sections of particles (i.e. each particle absorbs/scatters the power cri) . 
Thus, the first term in the right hand side of Eq. (2.1) reflects the decrease 
in I due to absorption and scattering by particles in volume ds. The second 
term represents the portion of the specific intensity incident on this volume 
from all other directions due to scattering from particles outside this volume. 
The third term represents the increase in I due to emission from within the 
volume ds. 

To calculate J one defines a so called phase function p (§,§’) which is 
the probability that radiance at ^'in a direction will be scattered into 
a solid angle about 's. It is defined by 




r) 


1 


( 2 . 2 ) 




Aw - 



where f is the scattering amplitude. J is then given by 


(2.3) 


J c ^ 


J- 

qn ) 
qn 


•^The name "phase function" has its origin in astronomy where it refers 
to lunar phases. It has no relation to the phase of a wave 


(2.4) 
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where the Integration over all w’ is taken to include the contributions 
from all directions s’ . 

In the above equations, the particle density and size can be at different 
locations and therefore p o anii p can be functions of r. It is sometimes 
convenient to measure the distance in terms of non-dimensional "optical” 
distance t defined by 

T c j 

With this definition, from Eqs. (2.1) and (2.4) 

at ^ 

-}- ^- 1 ^' V (2.5) 

Tliis is the basic integro-differential equation which needs to be solved 
for I for a vegetation canopy. The solution involves two major steps; 

(1) Calculation or Specification of the phase function p(s,s^ in terms 
of vegetation canopy properties. For any applications, this is a rather 
difficult task. It is even more difficult for the canopy because vegetation 
is extremely heterogenous and complex and the canopy can not be treated either 
as a regular or completely random medium. Also, the scattering and absorbing 
elements of vegetation canopy, namely, leaves, stalks, flowers, etc. are very 
large compared to the molecules and the aerosol components of air, and are 
characterized by a relatively high absorption coefficient; about 0.85 for the 
photosynthetically active radiation (PAR) and about 0.15 for the near IR 


radiation. 
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(2) Solution of the Equation for a given phase function and boundary 
condition. One of the procedures for solving an Integro-dlfferentlal equation 
j.ike (2.5) Is to substitute an initial guess for I In the right hand side 
and then Integrate the equation, subject to the boundary conditions on 1^ 
to get a new I which is then used again in the i'ight hand side to get a new 
solution. This iterative procedure Is continued till the value of I does not 
change (within a desired accuracy) . For vegetation (unlike the top of the 
atmosphere) , the upper surface of the plant cover is exposed both to the direct 
specular radiation and the diffuse flux of the scattered radiation from the 
sky, leading to a somewhat difficult boundary condition. Also, the optical 
thickness of the plant v^iiiopy is substantially higher (8-10) than that for 
the atmosphere (0.2-0. 6) leading to a slow convergence of the iterative pro- 
cedure. 

One simplification which vegetation canopy provides is that the emission 
from within the canopy is negligible and hence ^ ( V.- )c <3 in Eq. (2.5). 

The unique features of the vegetation canopy have been recognized by 
the various investigators. Generally, simpler problems are solved by imposing 
abstractions on the shape or boundary of the medium and on the form of the 
phase ftinctlon. We will attempt to delineate them in the next section. Before 
we do so, we will like to note the following two points: 

(a) Canopy as a parallel plane medium: 

One of the simplifications made by the majority of crop canopy reflectance 
models is to approximate 'tedium (canopy) with a parallel-plane infinitely 
extended medium, i.e. the one in which the medium can be split up into distinct 
layers (one or more) in which the optical and structural properties are constant. 

In this case, the specific intensity is a function of one dimension z perpendicular 

* See however Ellenson and Amundson (1981) who found delayed light emission from 
soybean leaves exposed to sulfur dioxide and used it as a means to detect 
plant stress. 
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to the layer and angles 0 and ^ defining the direction of incident beam. 
For this case, ?t *.2 radiative transfer equation may be re-expressed as 


' ^ X 

c\'C 

1 /iiT 


where K 


( 2 . 6 ) 


(2.7) 


o 


U r Ci-S ly ax. - 5 ^ a z 

and T is the optical distance in the z direction ( 
Eq. (2.7) can be formally integrated to give 


X. 

r 


T 

'*•0 


( 2 . 8 ) 


Basically, this equation states that the upward (downward)' radiance at optical 
path T, is a result of the upward (downward) attenuated radiance at tq plus 
that scattered into the beam along the path between t and tq. 

Even for the simple geometry, no closed form solution has been found for 
a general arbitrary phase function and one has to resort to computer based 
numerical solutions. Various canopy reflectance models either make some 
further approximations or find numerical solutions. 

(b) Polarization Effects ; 

The above formulation of the radiative transfer theory excludes the 


polarization effects. Such effects may be important for some classes of 
materials, e.g., vegetation canopies with waxy leaves - pine needles, rhododendron. 
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holly - which often produce a strong specular reflection or glaze. 

To include the polarization effect one replaces I(r,s) by a vector 
I(r,s,t) whose components are the so called Stokes parameters 1^^, 1^, U and V 
(Ishlmaru, 1978a) defined by 


■ h ■ 

(2.9a) 

2Re (Ej^ E 2 ) “ 2 ^osf 

(2.9b) 

2Im (E^ ^2 ^ “ ^^1^2 ^ 

(2.9c) 


where 

^ ^ ’ ^2 “ ^2 ^ ^ (2.9d) 

are the phasor representations of the electric field components and E^ 
in the mutually perpendicular directions (E^ = a^^ cos (wt - kz + a^) , 

Ey = a 2 cos (wt - kz + ^2 )) ^ plane wave propagating in the z direction. 

Whether one should stress polarization effects or not in the modeling 
effort is not clear at this time. Egan (1968, 1970) suggests that discrimi- 
nation potentials may exist in the asymmetric depolarization effects as a 
function of view angle. Curran (1981) argues that the polarized visible 
light is an indication of scene roughness which is a function of vegetation 
amount. He also finds a linear relationship between polarized visible light 
and vegetation amount. 

Vanderbilt (1980) had modeled plant canopy polarization response and 
Vanderbilt, Biehl, Robinson, Bauer and Vanderbilt (1980) have measured the 
linear polarization of light by wheat canopies. Tliesc scant but difinitive 
works do suggest investigating polarization effect&t especially in the context 
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of ground based reflectance measurements. 

We now move on to the discussion of the various crop reflectance 


models. 
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III. ANALYTICAL CROP REH^ECTAKCE MODELS 

We will now provide an overview of various crop reflectance models, 
analytical ones in this section and the numerical ones in the next section, 
with emphasis on the key assumptions and the agronomic variables used in 
the models. In Section V, we will compare and contrast various models and 
recommend a strategy for future work. 

(1) Models Dirt ctly Based on Kubelka-Munk (KM) Theory 

As an alternative to the numerical solution of the transport equation, 
in 1931, Kubelka and Munk proposed an approximate theory for a parallel plane 
random medium (See (a) of Section II) . The light flux is described in terms 
of two nonochromatic fluxes E_ and E_j_, travelling in the downward and upward 
directions, perpendicular to the plane of the medium. The variations of these 
•fluxes within the medium are described by means of two parameters a and y 
which are absorption and scattering coefficients, respectively, assumed same 
for E_ and E^ . More specifically, 

dE_ (2.1.1a) 

~ =-(a + Y)E+YE, 

d(-r) 


= - (a + y) \ + Y E_ (2.1.1b) 

The first term in the right hand side of Eq. (2.1.1a) states that the down- 
ward flux decreases due to its absorption and scattering, wihile it increases due 
to scattering of upward flux. Eq. (2.1.1b) has similiar physical interpretations. 

The relationships between fluxes E and E_^ and specific intensity I 
introduced in the preceding section are 
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n 

^ j J S>»0 (2.1.2b) 

The attractive feature of the KM theory is that it requires no more 

than simple algebraic operations to calculate the flux at any point in the 

medium. This attraction is in a large part responsible for the application 

of the KM theory in many applications of the radiative transfer theory. 

The drawbacks of the theory are the requirements of empirical determination 

of coefficients a and Y > and lack of theoretical basis and precise understanding 

of range of its validity. (See also Kortum, 1969). 

Subsequent to the initial proposal, a number of authors refined the KM 

theory and compared it with experimental data. The general conclusion of 

these investigations is that this two flux (E and E ) theory adequately 

“ + 

describes the experimental data if 

. illumination is diffused (isotropic light) and 

. medium is dull so that the light is diffusely scattered (ideal diffusing 
medium) 

It is not applicable if 

. a collimated or specular beam is incident upon the medium 
. medium consists of several distinct components that result in anisotropic 
reflectance. 

In connection with vegetation canopy, it should be noted that for a cloudy 
day, the canopy is illuminated by sun light which is predominantly diffused. On 
a clear day, the sun light is predominantly collimated. However, it becomes 
diffused very quickly as it enters the canopy especially when the canopy is 
reasonably dense. Further, the canopy does have distinct components which lead 
* E_ and E_j_ at the top of the plane parallel meditun are also knoxm as irradiance 
and radiant emittance (or radiant excitance) , respectively. 
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to a Isotropic reflectiince e.g., the reflectance difference between upper 
and lower surfaces of many plant leaves, bidirectional scattering effects 
of Individual leaves, inhomogenous distribution of leaf orientation, etc. 

The case o£ collimated beam can be included by adding two more fluxes, 
upward (F^) and downward (F_) collimated fli*:.js, resulting into the so called 
KM four flux theory (Ishimaru, 1978a). The equations which describe this 
extension are 

« - (t: + Y) E_ + YE^ + F_ H- F^ (2.1.3a) 

» - (a + y) E^ + yE_ + \ + ^2 (2.1.3b) 

» - (k + ) F_ (2.1.3c) 


«* - (k + + S 2 ) F^ (2.1.3d) 

where k is the absorption coefficient for collimated beams, and and $2 are 

the cocffioleuts of scattering from a collimated beam into a diffuse flux in 
the same direction and in the opposite direction, respectively. 

If one assumes that the collimated flux is only in the doraward direction 
(F^ = 0) the four flux KM theory is a three flux theory and is called as the 

Duntley theory (Duntley, 1942) involving 5 parameters. 

The Inclusion of an isotro pi c canopy reflectance in the theory is rather 


dE 


d(-T) 


dE. 


d(r) 


dF 


d(-T) 


dF. 


dr 
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difficult. One a'spect, namely, asymmetry between optical properties of the 
two sides of the canopy components, can be included by using different absorp- 
tion and scattering coefficients for fluxes in the upward and downward directions. 
For example. Park and Peering (1982) (see also Bunnik and Verhoef (1974)), for 
the case of diffuse flux only, modify the KM equations to the following equations: 

dE_ 

~ = (a_+ Y_) E_ + 

d( -T ) 

dE 

= - (“^+ Y+ ) E^ + Y_ E 

d( -T ) 

where a_ and are absorption coefficients associated with E and E^ 

respectively and y_ and the respective scattering coefficients. 

Before we discuss the applications of the KM theory to vegetative caiiopy 
reflectance, we should note the following: 

. (a) The hemi-spherical reflectance of the canopy can be calculated by 
forming the ratio of the upward directed diffuse flux to the downward directed 
flux (specular + diffused ) at the top of the canopy. Here it should be noted 
that since the diffuse flow within the canopy is assumed to be isotropic, the 
canopy reflectance is presumed to be Lambertian. 

* A Lambertian surface is one for which the specific intensity I (r, s) is 
independent of the direction s, i.e. radiation is isotropic. For this case 
from Eq. (2.1.2a) 

E_^ = I (t) 2ir l ^os^e 

i.e. the ratio' of irradlance to specific intensity (also known as spectral 


0 = 0 


= irl (t) 


0 =ir/2 


(2.1.5) 


(2.1.4a) 


(2.1.4b) 
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(b) Variable t can represent a canopy variable other than optical thickness 
(provided this new variable either is or assumed to be proportional to 
optical thickness) with accompanying slightly different interpretation of 
parameters. For example, Park and Peering (1982) set x as biomass per unit 
area. Allen, Gayle, and Richardson (1970) set x equal to cumulative leaf 
area index (LAI) defined as cumulative one-sided leaf area per unit ground 
area from the top surface of a stand to a plane at a given distance above 
the ground. 

Coming back to the specific use of the KM theory in connection with 
vegetative cano py reflectance, it seems that investigators at the U.S. 

Department of Agriculture at Weslaco, Texas were the first ones to make 

such use of the theory. Allen and Brown (1965) applied it to radiation in 

corn canopy and noted its inadequacy to account for the specular component 

of radiation and for the observed variation of plant canopy reflectance 

with the sun angle. A good test of the applicability of the KM theory 

was carried out by Allen and Richardson (1968) . They measured the reflectance 

(R) and transmittance (T) of stacked mature cotn-on leaves (normal aixd dehydrated) 

over the spectral range 0.5 — 2.5p and found that the deviation between theory 

and experiment is only about 1% , assuming certain wavelength (X) dependence 

of a and y (i>^ their notation K and S) for cotton leaves. They also gave an 

explicit solution to the KM equations for an actual plant canopy. 

In order to account for illumination of canopy by direct sunlight and 
phenomena produced by sun angle, leaf orientation, or other attributes, 2 
parameters KM theory was replaced by 5 parameters Duntley theory (Allen, Gayle 
and Richardson, 1970). Recalling that this theory assumes no upward specular 
light, its application assumes that the specular light i:\cident on the leaves 
as well as on the soil background is reflected as diffuse light. The effect 
of the sun zenith angle, 0, is included by assuming a sec 0 dependence of 
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specular light attenuation coefficient, k + + S 2 (see Eq. 2.1,3c). The 

Duntley equations fit the near infra-red experimental data on com canopy 
very well. 

In Table 2.1.1 are given the values of various fitted parameters that 
specify near-IR irradiance in a corn canopy. 

Table 2.1.1 

Fitted parameters that specify near-IR irradiance in a com canopy 
(Allen, Gayle, and Richardson, 1970). 


Parameters 

Fit 1 

Fit 2 

a 

0.000 

0.035* 

Y 

1.369 

•k 

0.736 

k 

0.000 

0.125 

^1 

0.609 

0.281 

^2 

0.978 

0.297 

Standard deviation 

3.2% 

3.7% 


* Laboratory values for single leaves at X = Ip. 

Fit 2 corresponds to the case when the laboratory measurements of optical 
constants, a and y, for a single corn leaf are used as constraints. Fit 1 
corresponds to no constraint case and suggests that the best fit to the near 
infra-red transmittance occurs when zero absorptance is assumed for the canopy 
(a = 0 , k = 0) . 

More recently, the KM two flux theory, with unequal absorption as well as 
scattering coefficients for upward and downward flux, has been applied by Park 
and Deering (1982) to the observed reflectance at 0.68p as a function of dry 
biomass for alfalfa and shortgrass prairie canopies (for alfalfa different 
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biomass levels were created through selective thinning within a larger uniform 
stand of alfalfa and the observations were made under various sky and illumination 
conditions). They tried to fit the data to several models reflecting various 
relationships between four parameters andy^ (eq. (2.1.3)) . For alfalfa, 

no noticeable differences were found in the fits for the following two models 

Model A: r = -/y^ >> 1 (Downward flux scattering much more than that 

of upward flux 

Mb.del B: a a a « a , y*y = y , (i.e. KM 2 flux theory) 

T* — — 

For shortgrass prairie no definitive best model could be established presumably 
because of the large scatter of the data. 

They also found that 

. optical parameters varied for differing illumination conditions and for 
sure depended on the sun angle. 

We should add that they found the general model (with four parameters a_, a^, 

y , y . ) "inefficient for computation" and hence did not use it in their analysis. 

T 

It should be emphasized that the best fit parameters for alfalfa and shortgrass 
prairie .canopies were quite distinct from each other (Sec Table 2.2.2). 

Table 7.2,2 


Parameter Alfalfa Shortgrass prairie 

K(10^ hectacre/kg ) 444-741 18.6 

R .019-. 033 0 

V 

(Oil reflectance, .165-. 258 =.186 

D («_-“+ + Y_ - Y+ )/2 

K ^ (D^ + (<_ + o<- V+ + f- )l/2 

+ Y_ + ) /(2y^. )|' " (K/y+) 


and hence 


. The KM theory based models, in priciple, could be used for automatic 
crop identification and possibly also assessment. 

We should point out that though the KM model has been explicitly studied 
when canopy has only one layer, i.e. it can be assum&d to be a. homogenous layer, 
it can easily be extended tor the case of many layers, with each layer being 
characterized by its own set of parameters. This extension requires ensuring 
the continuity of upward and downward fluxes at the boundaries of various 
layers as has been done by Chance and LeMaster (1977) in connection with the 
study of the Suits Model which we will describe next. 
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2. Suits Model and Ita Refinements 

In the KM theory, Duntley theory or In spirit equivalent model due to 
Allen, Gayje and Richardson, since the upward diffuse flux Is Isotropic 
or Lambertian by definition, the calculated reflectance Is Independent of the 
view angle while It does depend on the solar Incidence angle. Also, canopy 
geometry Is not taken Into consideration. 

The first canopy reflectance model that Is truly bidirectional, l.e. 
dependent on both the solar zenith angle and on the zenith angle and on the 
zenith angle of observation was developed by Suits (1972 a,b,c; Suits and 
Saflr, 1972). This was achieved by solving radiative transfer equation more 
exactly yeilding more realistic non-Lamb ertlan canopy radiance. The model 
also expressed parameters of the KM theory in terms of parameters, defining 
the canopy geometry and optical properties of the canopy components, which 
can be measured in the field. Canopy layers containing mixed components (e.g., 
leaves, flowers, and stalks) can also be included. In Suits' model, the canopy 
is idealized as a hrmogenous mixture of horizontal and vertical diffusely reflec- 
ting and transmitting flat panels that are considered to replace the original 
leaves by taking their horizontal and (two orthogonal) vertical projections. 

For this simple geometry, the KM model parameters are expressed in terms of 

* 

the following parameters 

Oj^ = average area of the projection of a leaf (or another vegetative component) 
on a horizontal planne. 

a = average area of the leaf wTien projected onto two orthogonal vertical 
planes, ( see Fig. 3.1) 

nj^ = number of horizontal projections per unit volume, 
n^ » number of vertical projections per unit volume, 
p = hemi-spherical reflectance of a leaf. 

T = heml-spherical transmittance of a leaf. 

0 = polar angle for incident specular flux. 

^Suits does not give the explicit derivation of these expressions in his publica- 
tions. See Slater (1980) for these derivations. 
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The horizontal projection o* and the two vertical projections o, of a leaf. 


Figure 3.1 
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p and T are assumed to be the same for both sides of a leaf. 

If more than one kind of vegetative components exist in a canopy layer, 

then the values of the KM model parameters are obtained for each type seperately 

and then added together to obtain the value for the layer. 

The radiative transfer equation is approximately solved by making an 

initial guess for the total radiance field I. For this purpose I is assumed 

as a sum of upward and downward diffused field and an attenuated downward 

specular fields. 

I * I , + I + I 
+ — s 

This initial guess is obtained by solving the resulting simplified radiative 
transfer equations which are the KM 3 flux or Duntley equations. The phase function 
is also assumed to be as a sura of three parts (u, v, w’; Suits 1981 ) corres- 
ponding to these three fields. Each of these phase function is explicitly 
expressed in terms of the canopy goemetric parameters and the optical proper- 
ties of the canopy given earlier in this subsection and the polar view direction 
0^ and azimuthal angle between sun and view direction . The initial guess 
together with the phase function is then used to calculate the source function 
K in Eq. (2.7). This value of K is substituted in Eq. (2.6) which is then 
solved to get an updated estimate of I along a particular direction. In principle, 
this iterative procedure could be used to generate a solution to any desired 
degree of accuracy. However, Suits stops at the first iteration. 

The procedure leads to a closed form expression for canopy reflectance leading 
to easily implemented computer codes. (There is now available a code in BASIC 
language which can be run on an Apple II microcomputer). Suits (1972a; see also 
Slater, 1980), and Chance and Cantu (1975) give explicit expressions for the 
canopy reflectance for single layer and two layer canopies. The only additional 
parameters these expressions have are the soil reflectance p , thickness of 

* The relation for w* had a missing factor in the original publication (Suits, 
1972a). 


- 21 - 


various layers in the canopies and re'lative fractions of diffuse and specular 
incident lights. Also, the parameters o^, nj^, n^ occur only in the 
combination 

The Suits model has been tested against experimental data by Suits and 
his colleagues at his institution (Suits, 1972a, b, o; Suits and Safir, 1972; 
Colwell, 1974) and by many other investigators at other institutions, most 
notably at Pan American University (Chance and Cantu, 1975; 

Chance, 1977; Chance and LeMaster, 1977; 1978; LeMaster and Chance, 1978; 
LeMaster, Chance and Wiegand, 1980), and at Netherlands Interdepartmental 
Working Community for the Application of Remote Sensing Techniques, NIWARS 
(Bunnlk & Verhoef, 1974; Verhoef & Bunnik, 1975, 1976, 1981; Bunnik, 1978). 
These studies provided a very good understanding of alterations in the 
reflectance caused by the botanical features and the geometry of the inter- 
action (Cause-effect relationship) . They also pointed out the shortcomings 
of the Suits model and refined and extended the model. 

It is beyond the scope of this report to go into the details of the 
various insights in the cause-effect relationships and on the shortcomings. 

The readers are referred to the works cited in the preceding paragraph, espe- 
cially works by Verhoef and Bunnik (1975) and by the Pan American University 
group. Here, we only point out that there is, in general, a good agreement 
between the model and the experimeiit.al data for spectral reflectance, and the 
effects of variations in leaf area index, average leaf inclinations and soil 
moisture on the reflectance could be simulated very well. For X = 600-690 nm, 
with small penetration, a single layers models fits well, while for X> 690 nm, 
with deeper penetration, a multilayer model is required. For X = 500-600 nm, 
with moderate penetration, a two layer model is the best choice. (Chance and 
LeMaster, 1977). Here, we should point out that Chance and Cantu (1975) 
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developed the mathematics which enabled the extension of the Suits model 
for an arbitrary number of layers. Also, Bunnik and Verhoef (1974) intro- 
duced unequal optical properties at the upper and lower sidch of the canopy 
material. 

The poorest agreement between the model and the experimental data is 
found for very early and very late in the growing season when the green plant 
biomass is low and the ground cover is Incomplete (LeMaster, Chance and 
Wiegand, 1980). At these times, the assumption of continuous and uniform canopy 
of the Suits model are not well satisfied. There is another situation where the 
plant canopy is not uniform and homogenous horizontally and the Suits model does 
not agree with the observation. This situation leads to the so called "row 
effect" and is described below (Suits, 1981). 

Many crops are t/lanted. in rows by machinery. Uuon emergence of the plants, 
the bare soil between rows is still the dominant feature which reflects incident 
daylight. As growth continues, the vegetation grows both higher and spread out 
over the inter-row regions covering the bare soil. For a considerable time during 
the early part of the growing oeason, the strips of bare soil between rows and the 
increasing density of vegetation along the rows become equally important in their 
contributions to canopy reflectance. During this time, the direction of sunlight 
relative to row direction will change the relative influence of vegetation, and 
bare soil. Ivfhen the sun is directed along the row direction, the bare soil 
is fully illuminated, but, when the sun is directed across the rows, the soil is 
largely in the shadow of the standing vegetation along the rows. Thus, Landstat 
can receive different alterations due only to the way the rows trend relative to 
sunlight. An inference that such layered radiation is due to a change in some 
important agronomic feature, could thus, be in error. 

Field measurements of wheat (Jackson et al, 1979a) and soybean (Vanderbilt 
et al, 1981, Randon et al, 1982) show that row direction relative to the sun light 
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doea ohmiRc llvt* dmractev of the acnttereil nml reflei'tod radiation. In fact, 
the rcflcctmtco may change by a factor of 2 or 3 for the aamo canopy dept^t'ding 
upon the direction of rowa relative to Hunlii^ht. 

Verhoof and Munnik (1^)76) were the first ones to extend the Suita model 
to Include the effecta of rows by nssuminjt rectanRulnr crosssectiona of crops 
planted in rows, with random avrauRement of leaves in each row and onlv soil 
between rows. AccordiUR to Smith and Ranson (l‘)70), they undertook a detailed 
Rcometrienl analysis of the canopy phase function relative to direct solar 
flux, with shading allowed. However, the approximation la made that the diffuse 
flux can be treated ns In the homogenous .Suits model. For both types of fluxes, 
an appropriate view probahllity function is developed that is consistent with 
the row atruetore. The soil contribution to the canopy reflectance la carefully 
developed consider Inc, the row structure, hvjt even then the contribution Is 
dlaconl Inons. Leaf angle distributions are handled via vertical and horizontal 
projections as in the basic Sult.s model. The row model was verified against 
spectral measurement on wheat and It did predict the angular dependence of ref- 
leetanee relative to viewing aslmuth. 

More recent Iv Suits (1981) Included the row effects In his model for mon> 
general cross-sections of crops by essentially multiplying the KM parameters 
In his model by a row modulation factor depicting the geometry of the rows. He 
applied the model to v^heat and found that the new model does give the result 
slmUlar to those of field measurement and to those obtained by Verhoef & lUmnlk 
(197b). 

LeMaster and Chanee (1978) fovmd another disagreement between the model 
and experimental data. For Feivlamo wheat, tlu' observed reflectance at 550 and 850 
nm Increased with the sun zenith angle, while the Suits model predicts a decrease. 
Slmlllarly, for a fixed sun angle, the reflectance Increased with Increase in the 
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Observer zenith angle, contrary to decrease predicted by the Suits' model. 

This disagreement may be due to the assumption of Lambertian reflecting 
leaves used in the Suits' model for specular irradlance (LeMaster and 
Chance, 1978) and/or the assumption of no oblique leaves (Verhoef and 
Bunnlk, 1981) . 

Suits model has been extended by the Dutch group (Verhoef and Bunnik, 

1981) , to include scattering and extinction functions for canopy layers 
containing fractions of oblique leaves according to a given leaf inclination 
distribution function (LIDF) . In this model, called SAIL (Scattering by 
Arbitrarily Inclined Leaves), the Suits model simplification of canopy geometry 
to exclusively horizontal and vertical panels is replaced by discretization of 
canopy LIDF to a set of frequencies at distinct leaf inclinations 0, . SAIL 
model seems to have removed the limitations of the Suits model and can predict 
crop reflectance for most of the leaf inclination distributions found in 
nature. 

Another extension of the Suits model was carried out by Sadowski and 
Malila (1977) to forest canopies. They introduced slope angle and additional 
parameters characterizing ^;round surface. Their model is empirically calibrated 
and is applicable, for a uniform canopy; it is not driven by tree size and 
spacing. The model lacks an explicit geometric component to account for shadowing. 
The model depicts fairly accurately the observed reflectance, including the phe- 
nomenon that open canopy stands under low angles of solar illumination will 
have reflectance similiar to a closed canopy stands under higher illumination 
angles (Strahler, 1981). 

We conclude the discussion of the Suits model by pointing out that Chance 
and LeMaster (1978) proposed a light absorption model (LAM) for vegetative plant 
canopies based on Suits canopy reflectance model. Both of these models have 
the same set of experimental parameters. Reflectance model has its value In 


in terras of characterizing the canopy, while once this has been done, the 
absorption model can be used to determine the absorption of light in the 
photosynthetically active region (PAR) of the spectrum (4000-7000A® ) 
which in turn is required for determining the crop growth. LAM's predictions 
are found to be in agreement with the experimental data for Penjamo wheat 
in the soft dough stage (98 days from the emergence with an LAI of 3.5). 
Also, the variation of percent canopy absorption as a function of LAI seems 
to agree with the measurements on Plainsman V wheat. 


IV 
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NUMERICAL CROP REFLECTANCE MODELS 


4 


In this section, we will briefly discuss those models in which numerical 
means are used to either solve completely and directly the radiative transfer 
equation (2.5) or (2.8) or part of it or its equivalent manifestation. The 
hope is that one will then be able to analyze more realistically the canopy 
reflection for a larger verlety of canopy geometries. At the outset, we 
should point out that these models and their intricacies are hard to grasp even 
by those who have developed competitive approaches. This together with the 
difficulty in obtaining the computer codes which implement, various models 
are in a large part responsible for almost an order of magnitude less use of 
these models as compared to the models described in the preceding section. 

(1) Smith and Oliver Model 

It is one of the more known numerical models proposed essentially at the 
same time as Suits model (Smith and Oliver, 1972; 1974; See also Slater, 1980; 
and Smith and Ransou, 1979). Here, a direct attack on the numerical solution of 
the radiative transfer equation (2.8) for a plane parallel medium is made. 

In this model, the flux within the canopy is allowed to propagate in 
discretized (0*,‘l>’) directions rather than only in upward and downward direc- 
tions as for the KM thoery based models and Suits model. 

The main feature of the model which allows for this generality is the 
calculation of the layer phase functions from the angular distribution of the 
foliage elements and the reflectance and transmission properties of these 
elements with respect to the discretized 0*, <J>^ source directions. A foliage 
element inclined at an arbitrary orientation with resv*^;ct to the source direction 
permits^ according to the Lambertian response, the scattering by transmission 
and reflection of the incident flux to upper and lower hemispherical sectors. 
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For each foliage Inclination represented in the canopy a set of integration 
limits on the scattered radiation from a scatterer is defined. For a given 
layer the distribution of flux is then weighted by the frequency distribution 
of foliage inclinations occurring within the layer. 

The vnodel initiates an iterative solution of Eq. (2.8) by using the 
zero order flux above the canopy to generate via the phase function of the 
first layer the estimated flux in layer one. This is then used together 
with the phase function for layer two to calculate the estimated flux in 
layer two and so forth. Subsequently, reflection from the soil boundary 
generates upward moving flux, again in a *set of discretized 0*, directions. 
Frocosslng ‘is continued until all flux levels within layers reach equilibrium 
values . 

The reported version of the model is a Monte Carlo or stochastic 
implementation of the above processes. The following description adapted 
from Smith and Oliver (1972, 1974) and Slater (1980) gives more insight into 
the model . 

The direct solar radiation is treated as a set of independent source 
vectors to the canopy. For simplicity, all the vectors are divided into 
source sectors formed by partitioning the herai-sphere into 10 degree 
inclination bands and further subdividing the bands to form 20 degrees 
azimuthal sectors. The midpoint of the sector is taken as the direction 
of the diffuse flux from that sector. The Interaction with the canopy of 
each of these Intial 181 radiation sources is treated independently. 
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The diffuse flux resulting from the interaction of global radiation with 
a canopy element or with the background become now sources which may further 
Interact with the canopy. The downward directed flux is combined with 
the appropriate hemispherical band of diffuse sky radiation. Upward 
directed flux is treated in a similar manner as diffuse sky radiation 
except the direction associated with each sector is the opposite from 
Incoming radiation from that sector. 


Geometrical measurements of a canopy provide a frequency distribution 
of foliage inclination angles. These can be integrated using Simpson's 
rule to obtain a cumulative frequency distribution that is nomallzed 
and partitioned into areas of equal probability. The interaction probabi- 
lities within each canopy layer are then calculated for both incident and 
emergent fluxes at each specified source direction. The following expression 
due to Idso and DeWit (1970) is used in the model: 


P 


0 



S X 


sin 0 j 


(4.1) 


PHIT - 1 - Pq , 


(4.2) 


where Pq is the probability of a gap, is the probability of an 

Interaction, 0PM is the mean canopy projection in the direction of the 
source, 0 is the source inclination angle, and S takes b value from 0 
to 1 depending on the density of the canopy. It is usually adjusted by 
optimizing the reflection prediction for either one view angle or a set 
of wave length since it does not change with either wave length or view 
angle. For large LAI, PHIT is Insensitive to S and S =* 0.1 is 
satisfactory. 
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A random number generator is used to generate a random number. A given 
source finds a gap in the top layer of the canopy if a random number is 
smaller than PHIT. The flux in this direction passes through the top 
layer and reaches the next canopy layer. The absence of a gap necessi- 
tates the determination of material type with which contact has been made. 
This is accomplished by sampling from the distribution of canopy cons- 
tituents. The orientation of the leaf is determined by sampling from 
the inclination distribution and a uniform azimuthal distribution. These 
parameters determine the direction cosines of the leaf from which 
the angle between the leaf and the source is determined. The optical 
properties of the leaf are then utilized to calculate the flux exiting 
the leaf in all directions. Each sector of the hemisphere on the reflecting 
side of the leaf receives flux for each wavelength according to the equation, 

I =1^ P sin (Qj^g) (sin^ " sin^ 0^^) (^*3) 

l8 

= source spectral flux 

p = material spectral reflectance 

= angle between the leaf and the source 

02 , Gj^ are the inclination angles defining the hemispherical 
band 

The solid angle sectors receiving reflected and transmitted flux are defined 
in the same manner as for canopy flux sources only extended to include the 
entire sphere about the leaf. The leaf is not necessarily horizontal so a 
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sector receiving reflected flux from the leaf is not necessarily directed 
upward with respect to the local vertical. Hence » the direction cosines 
of the flux sector are rotated to the local vertical system and the flux 
pooled with the flux in the appropriate source band. Transmitted flux 
is calculated and treated in the same manner as reflected flux. 

Flux which passes through a gap or is reflected or transmitted downward 
from an upper layer of the canopy interacts with the next lower layer. 
Flux which reaches the soil sureface is reflected into each of the upward 
directed source bands according to the equation: 


where 


I = I 


sin (0) 


(sin^ 0, 


- sin 




( 4 .«) 


= soil spectral reflectance 
0 = source inclination angle 


Upward directed flux from a lower layer of the canopy or from the background 
reaches the next higher layer and may interact with it. The upward directed 
flux from the top layer escapes the canopy. 


The interaction procedure continues until the level of flux in any source 
direction within any layer is below a threshold value. The flux exiting the 
canopy into each of the bands is separately accumulated. The ratio of the 
flux intercepted by a sensor placed within one of these bands to the global 
radiation intercepted by a vertical sensor with the same filed of view gives 

I 

the canopy directional reflectance. 


Several iterations through the model will reduce the standard error of the 
predicted reflectance. For each pass, it may be desired to generate random 
vectors, e.g., constituent optical properties or global radiation flux, to 
drive the model. The model is capable of generating random vectors from a 


multivariated normal distribution or will use a mean vector on each pass 
if desired. 

Detailed calculations were made for the case when the canopy is 
modeled as a three-layered vegetatior ensemble, with each layer containing 
two types of vegetation elements which are assumed to be Lambertian 
scatterers. The model was evaluated (Smith and Oliver, 1974) against 
measured data for shortgrass prairie vegetation with a leaf area index 
of 6.5. For a vertical view angle, agreement is good except for the 
chlorophyll absorption band x^here the model predicts lower reflectance. 

Smith and Oliver suggest that this discrepancy may be due to their assump- 
tion of the surfaces being Lambertian whereas there is some indication that 
the reflection characteristics are non-Lambertian in regions of high absorp- 
tion. Off-angle predictions of the model are qualitatively correct, but do 
not display the same precision as the vertical view case. The comparison 
also suggests that the Smith and Oliver model compares more favorably with 
the experimental result than those based on the Kubelka-Munk theory. 

The model has also been used to understand the observed wide variability 
in diurnal reflectance (solar zenith angle dependence) for lodgepole pine 
and two grass canopies (Kimes, Smith, and Ranson, 1980). The model correctly 
simulated this variability and suggests that this variation is caused by 
variations in anisotropic sky irradiance, canopy component geometry and optical 
properties, and the type of reflectance measurement. 

A major difficulty of the current implementation is the pooling of the 
outgoing radiance into 9 directions only. That is, outgoing azimuthal directions 
are averaged. It should be noted, however, that incoming source azimuth direc- 
tions are included. The model requires considerable time for the Monte Carlo 
analysis of a canopy. In principle, however, the approach is an accurate 
representation of the radiative transfer processes occurring within a plane- 
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parallel medium. A deterministic version of the model that Included 
outgoing azimuthal dependence could greatly enhance the utility of 
this approach. 

(2) Ross and Nilson Model 

This model is again for plane parallel canopy and has been deve- 
loped by two Estonian investigators Ross and Nilson(1966; see also 
Ross, 1975 for other references). They consider plant canopy as a turbid 
anisotropic and statxstlcally homogenous (in the horizontal direction) 
layer and the radiation field inside the canopy is determined by numerically 
solving a modified radiative transfer equation. 

They describe the geometrical structure of the canopy as a collection 
of numerous thin horizontal laminae. The thickness of each lamina . is 
selected so that a ray propagating in the vertical or nearly vertical direc- 
tion intercepts only one vegetation element; in this way the number of 
scattering events in each laminae is at most 1. 

They further assume that each lamina has sufficiently large horizontal 
surface such that it contains a large number of vegetation elements and one 
can assume a statistical distribution of the spatial orientation of the leaves 
over that surface. This distribution is assumed to be the same for all 
laminae. 

The vegetation canopy geometry is characterized by: 

(a) Uj^ (z) : the amount of vegetation elements of type k in unit layer volume, 
at distance z in the canopy. The leaf area index Lj^ for vegetation element 
of type k is related to Uj^ by 
h 

I 

0 


Uj^ (z) d z 


(4.4) 


where h Is the height of the canopy. 

(b) gj^ (z,T')5 the distribution density of the normal to the top surface 
of the vegetation elements of type k In a direction r. This Is used to 
calculate the so called G function which may be Interpreted as the mean 
projection of a unit foliage area In a particular dlrectlonS^ yhey derive 
a radiative transfer equation not for a volume element but for each 
elementary lamina using Intensities averaged over the coordinates x and y 
( z Is the direction perpendicular to the canopy) . This Is the equation 
which Is solved to calculate canopy reflectance. 

In their formalism, the optical properties of the vegetation elements 
are defined by general wavelength dependent scattering functions. However, 
because of the difficulty In solving the general problem, In detailed 
calculations, they assume identical scattering properties for the two leaf 
surfaces and use reflection and transmission coefficients as the optical 
property characterizing the vegetation elements. They also assume that 
gj^ Is Independent of z, l.e. distribution function of the leaf orientations 
Is almost Independent of z. 

They Investigated In detail three cases of leaf orientations J(l) all 
leaves horizontal, (2) all leaves vertical and (3) leaf normals uniformly 
distributed. For case (1) the radiative transfer equation can be analytically 
solved with reflectance exponentially decaying T/lth LAI. For other cases, 
the equation has to be solved numerically. 

On the basis of limited available published work in the english language, 
it appears that the model seems to agree well with the experimental results. 
However, because of the assumptions on basic laminae, the model can not take 
into account row effects without some drastic modifications. 

(3) Adding Method 

Very recently Cooper, Smith, and Pitts (1982) proposed a model based on 
a so called adding method (Van de Hulst, 1980), in which they do not directly 
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solve the radiative transfer equation. Instead, they divide the canopy 
into a set of appropriate sublayers, with each sublayer's reflectance 
properties characterized by two matrix operators R and R* (one for the 
top and the other for the bottom of the sublayer), and transmittance 
properties by operators T and T' . The elements of these operators give 
the flux transferred from one direction to another. Elements of R for 
the top sublayer will by definition give the canopy reflectance. 

Explicit expressions for these optical operators are calculated 
numerically in terms of the optical properties of the canopy components 
and geometrical properties of the canopy is assumed to comprise of indi- 
vidual Lambertian components characterized by reflection, transmission, 
and absorption coefficients. The components are assumed to be oriented 
isotropically with respect to the azimuth but a leaf slope distribution 
f(0j^) is specified. The other parameters used are LAI and reflectance 
of soil. 

The leaf angle distribution and LAI together with the Lambertian 
assumption, are used to calculate the proportion of the incident flux 
which does not interact with the canopy layer, probability of gap (see 
Eq. (4.1)) and to compute the total canopy projected area in the directions 
of incident and exitant flux. These calculations coupled with optical 
properties of the vegetation elements are used to calculate optical matrix 
operators . 

It should be noted that the size of various matrix operator and computer 
storage requirement will increase rapidly as the number of discrete ranges 
in which the incident and view directions are discretized is increased, 
(discretizing interval is decreased). In the initial implementation, the 
sun zenith and view angles ranged from 5 to 85 degrees and were discretized 
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in 10 degrees increment leading to 9 x 9 matrix operators. 

The model has been checked against the reflectance measurement for 
a Blue grama and soybean canopies and against the Suits model. It appears 
that the model generally gives good results at infrared wavelengths and 
overestimates the reflectance at visible wavelength for soybean, particularly 
as compared to the Suits model. For Blue grama it yields excellent results. 

The authors suspect that the discrepancy is because the model is sensitive 
to asymmetrical scattering properties of the leaves; in the simulations, 
p and T for the leaves were assumed independent of leaf side which is 
probably true for Blue grama but not for soybean. Also, for soybean the 
assumed form of probability of gap may be less valid than that for Blue 
grama . 

We conclude the discussion of this model by pointing out that the model 
does not take into account the row effects and at this time it is not clear 
how it could be Incorporated in the model. 

(4) Weinman and Guetter Model 

In this model proposed at the same time as the Suits model (Weinman and 
Guetter ,197'iJ) the canopy is also idealized as a plane parallel medium. The 
radiative transfer equation is solved using the method of discrete coordinates 
(Chandrasekhar, 1950; Lenoble, 1956) and four-point quadrature. The phase 
function is represented essentially as a series of Lfsgendre polynomials (series 
is truncated to include only the first four terms) . 

In their model they emphasized the need for the reduction of the number 
of parameters to characterize the optical properties of the vegetation canopy 
and the desirability of Including atmosphere in the same model while describing 
the penetration of solar irradiance. Their model requires only three parameters 
to describe either the canopy (modeled as a single layer) or the atmosphere. 
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These parameters are; 

a; albedo for single scattering, (ratio of scattering to scattering and 

absorption) . 

x: optical depth or LAI 

g: asymmetry factor which is a measure of the deviation of scattering 

from an isotropic scattering. 

They found that the data on com canopy reflectance at Ip wavelength 
used by Allen, Gayle, and Richardson (1970) can be fitted very well with 
their model with the following parameters: 


Layer 

height 

(m) 

T 

a 

g 

atmosphere 

2.5<z< 3000 

0.087>t>0.0 

l.O 

0.75 

corn canopy 

0.0<z< 2. ’5 

3.17>t>.08 

0.98 

-0.44 


(Note that the value of a “ 1 implies no absorption, which is consistent 
with the conclusion of Allen, Gayle and Richardson (See Section III). 

They emphasize that the three parameters can be derived from laboratory 
measurements of the optical properties of the leaves and canopy geometry. 

They also give the following relations between these parameters and those 
occuring in Duntley equations (eq. (2.1.3). 

2(l-a) =a/ (k+Sj^+S 2 ) (4.5a) 

I (1 - ag) - (1 - a) = Y/(k + + S^) (4.5b) 

dx = (k + Sj^ + S 2 ) d(LAI) (4.5c) 

This reviewer is not aware of any other comparison of the prediction of the 
model with other measurements. Therefore, its potentials and limitations 
can not be assessed. It appears that its extension to include the row effect 
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tnay be difficult, and may have to be done by sacrificing its simplicity. 

(5) Welles and Norman Model 

All the models discussed so far idealized the vegetation canopy as 
a parallel plane medium consisting of layers. We will now very briefly 
discuss a model which is truly three dimensional. It was proposed in 
1979 (Welles and Norman, 1979) and seems to be very promising. Unfortu- 
nately details of the model have not yet been published. Therefore, one 
is unable to make a definitive assessment of the potentials and limitations 
of the model. 

In this model, a canopy is modeled to consist of a finite number of 
a three dimensional geometrical figure like an allipsold or a cylinder 
which could be spaced in onu of many patterns - regularly spaced, densi''ly 
grouped (including overlapping) at regular spacing, sparse groups widely 
separated - to more realistically reflect a particular vegetation including 
those planted in rows. Within each of these geometrical figures, the foliage 
could be chosen to be distributed randomly or in a non-random fashion 
(e.g. different foliage density in the interior of the figure than on its 
periphery). This density distribution is chosen to possess the given foliage 
angle distribution. Once the foliage has been distributed in each geometrical 
figure, the attenuation of incident beam as it travers the collection of 
figures is calculated. This in turn is used to numerically calculate the 
diffuse radiation by using a scattering phase function, reflecting the 
optical properties of the canopy elements. In visible region, the assumption 
of single scattering is used while in infrared region- multiple scattering 
is allowed. 

The. next step which is rathe.r Innovative is to "transform" each point in 
the finite array of geometrical figures to an equivalent Infinite plane parallel 


medium canopy by choosing a depth In the plane parallel canopy that 
has the same diffuse penetration probability considering both the 
upward and downward radiation flujc. )'hat is, for an arbitrary point F 
in the array, a point H, in an infinite one layer canopy, is chosen, 
defined by its depth in the layer, such that the transmission at F in 
the upper (lower) hemisphere is the same, as upward (downward) transmission 
of flux at point H. Slmlliar equivalence is done for reflectance. Once 
equivalence has been established for all points in the array, the plane 
parallel medium is used to calculate bidirectional reflectance. SunHt 
and shaded leaves are treated separately. 

(6) Other Multl-dlmenslonal Models 

As plants grow, they possess characteristic shapes that govern 
the spatial arrangement of their reflective matter. Further, tnese shapes 
obscure varying amounts of soil or understory vegetation and also cast 
shadows on the soil, understory, or other plants. There are a set of other multi- 
dimensional models (as contrasted to plane parallel canopy case) in which 
goemetric form factors for the rows and an analysis of shadowing play a 
dominant part. These models employ primarily geometric optics and multiple 
scattering is ignored. Egbert (1976, 1977) developed a model which allows 
calculation of optical bidirectional reflectance from shadowing parameters 
of surface projections or perturbations. Strahler and Li (1981, Li, 1981) 
model low density timber stands as a collection of randomly spaced cones. 

Each cone has a fixed base/height rp' io, and is taken to be a flat Lambertian 
reflector which absorbs visible wave length differentially. Tree heights 
are assumed to be log normally distributed and tree counts from pixel to 
pixel vary according to Poisson distribution. This model can be inverted 
to yield tree heights and spacing from a remotely sensed reflectance data. 
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Other geometrical models have been proposed by Jackson et al (1979b) and 
Richardson et al (1975). 

These models have their role in the crop canopy reflectance modeling 


but perhaps not as much as the ones which include multiple scattering. 
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V. Sunmiary and Strategic Recommendations 

In the last three sections we have provided an overview of the 
radiative transfer theory and how it has been used to calculate the 
canopy reflectance given the properties of the canopy constituents 
and canopy geometry. Each approach to its use has led to a model for 
canopy reflectance. All the evidence collected, to support various 
models, suggest that radiative transfer theory, which is a macroscopic 
approximation to the interaction of radiant energy with matter, is appli- 
cable to crop canopy reflectance. Thus* it should continue to play its 
pivotal role in crop canopy reflectance modeling. 

Most of the models, with the exception of those based on the KM theory 
and Suits model, have been tested mostly by their authors using limited 
experimental data base. They do generally agree with the observations used 
in their testing (otherwise it is reasonable to conjecture that the model 
would not have seen the publishers ink) . To determine their range of 
applicability and relative merits and limitations, it is desirable that a 
uniform and extensive data base for 3 or 4 crops, in different stages of 
development and planted in different geometric configurations, be made avail- 
able for exhaustive testing of various models either by the authors or by 
a centralized testing group. The data base should include measured values 
of the agronomic, optical, and geometric variables characterizing the 
canopy as well as measured reflectance for a set of incident radiation angles, 
view angles, and atmospheric conditions. This testing if implemented will 
clearly establish which model could or should not be usdd in a given situation. 

Even without this test, it appears that none of the models will work 
for all crops and for all conditions. Therefore, parallel to the testing, 
development of crop reflectance models should continue. Here, though, in 
principle, it is desirable to develop as many models as possible, in practice 
one may have to choose only a few. If so, it appears that the best strategy 


will be to develop two types of models. 

One model should be simple and analytical (or at least semi-analytlcal) 
which Is easily comprehensible and could be used with limited computing power 
and time to calculate reflectance from known canopy properties for a fairly 
large set of crops and canopy geometries. A model adapted from Suits model 
seems to be a good candidate. This modified model should allow inclusion 
of unequal optical properties of vegetation elements on its txro sides, various 
leaf angle distribution functions, and row effects. As discussed in Section III, 
these inclusions have already been done seperately, but not in one encompassing 
and tested model. 

The other model should be a comprehensive one and be capable of incor- 
porating detailed properties of the canopy, without making any significant 
approjcimations in the process. This model should be characterized not by 
its simplicity and comprehendibility but by its accuracy in calculating the 
crop reflectance from canopy parameters. Such a model will, of necessity, 
require numerical solution of the radiative transfer equation and hence may 
not he very kind to computer storage and time requirements. The more difficult 
part of the development effort of this model is the characterization of 
scattering phase fi.vj<ction in terms of canopy variables and not the procedure 
for iterative numerical solution of the radiative transfer equation. The latter 
technology has been fairly well developed in many other applications of the 
radiative transfer equation. 

Since the simpler model is an adaption of several existing models, it could 
become available in a comparatively short term ( 1-2 years). The comprehensive 
model involves some original development and my require a longer term effort. 

All the canopy reflectance models proposed to date have a major deficiency. 
They do not include time as an implicit variable. That is, if one wants to 
calculate crop reflectance at different times (stages of development) , one has 
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to input in the model canopy variables for each time. Getting information 
for this input is rather time consumraing. Therefore, it is desirable 
that a reasonable effort be made to combine a canopy reflectance model 
(preferably a simple one like the Suits model) with a vegetation growth 
model. Such a model will then provide a natural link to the temporal profile 
models (Badhwar, 1980; Badhwar, Austin and Carnes, 1982; Badhwar and 
Henderson, 1981) which have been used for automatic crop classification 
and crop emergence date and growth stage determinations. 

To date most of the models have been used as a research tool and for 
understanding, i.e. for defining the proper Instrumentations (e.g. spectral 
bands of sensors), in interpreting data, for assisting in the identification 
of appropriate transform of reflectance in various wavelengths which may 
be insensitive to some canopy parameters, and for identifying potential 
causes of abnormal observations. 

They have also been the potentials for 'forecasting' reflectance 
for hitherto not tried set of canopy parameters. However, in light of the 
overall goal of crop identification and crop growth stage and quality deter- 
mination from the reflectance data, it is imperative that various crop reflec- 
tance models be investigated to assess their capabilities in correctly and 
uniquely determining the canopy parameters of importance like LAI, solar 
radiation intercepted (SRI), etc. from the reflectance data. In other words, 
they should be tested for their invertability. As part of this testing one 
should also determine how sensitive the results are to the variations in the 
reflectance data. Also, the invertability of reflectance for a given X vs. 
that for a linear combination of reflectance for many wavelength bands (e.g. 
Kauth-Thomas Greenness) should be investigated. Like with any such endeavor, 
first simpler models should be tested for inversion, including those which can 
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be semi-analytically inverted, followed by complex models requiring numerical 
inversion. Mathematical techniques for the inversion or what is more 
commonly referred to as systems identification are generally available but 
what is required is their adaption in the present context. 

In connection with the inversion, we should point out that it is possible 
that the present models, especially the more complex ones, may not be attractive 
from the inversion point of view. It is therefore, desirable to look into some 
reflectance models which are invertible but have not yet been tried in the context 
of crop canopies. One such nodel is the so called multi-thin layer model which 
has been used in designing optical system which modifies spectral composition 
of an incident radiation in a given manner (see e.g. Dobrowolski, 1981). 

As the name implies, in this model, the optical system is considered as 
a set of thin film layers system. The main physical properties that can be 
modified by such a system are transmittance, reflectance, absorption, and 
polarization. The system is designed to meet a required performance for one 
or more of these optical properties at a selected wavelength or in a certain 
wave length region. That is, the method can be used to determine the construction 
parameters of a thin film device, namely, refractive index n, absorption coeffi- 
cient k, and thickness d. The method has been demonstrated to be a practical 
one for most coating systems. 

To adapt the technique lor vegetation, for a given reflectance R(X, 9, (j)) 
which depends on wave length X, and view direction (6, (})) , one will use the 
method to calculate the parameters n(X),k(X) and d(X). One hopes that the values 
of these parameters will be sufficiently different for different crops to allow 
crop identification. It is conceivable that one may need to use R at different 
times and/or a multilayer system with a set of parameters n, k, and d which vary 
continousl> as a function of distance perpendicular to the thin layer. 
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At the outset It should be pointed out that there is aii important 
difference between the thin film devices for which the method has been 
developed and the canopy system. This may lead one to question its 
validity. A canopy structure causes an incident radiation to scatter in 
all possible directions instead of being reflected in the specular direction 
as for a thin film system. Ignoring this difference one may go ahead and 
obtain an equivalent multilayer system which can generate for a given Incident 
radiation the specified amount of scattered radiation in a given direction. 

For a different view angle, another set of parameters for the equivalent 
system may be obtained. If the two sets of parameters are close enough, 
then there is a high probability that as long as two crops have two different 
spectral characteristics, two different equivalent systems may be generated 
to tell them apart ( a successful crop identification technique). If not, 
one may have to modify the thin film technique to include scattering. 

We conclude this section by enumerating a set of modeling areas which 
are relevant to the crop canopy reflectance, but were outside the scope of 
this review. They should be looked into more d tail to determine a desirable 
strategy for further work. These areas are as follows: 

(1 Modeling of reflectance from single vegetation component, e.g. leaf as a 
function of wave length as well stress condition. This has the potential 
of allowing stress condition identification from the reflectance data. 

(2) Modeling of reflectance in the thermal infra-red region (3 -20 ym) . This 
is rather useful in light of measurements made in this band by Thematic 
Mapper on board Landsat-D. Also, it could lead to the assessment of water 
status of the bare soil as well as vegetation canopy form the reflectance 
data"whlch,if known during the early stages of growth, could be useful 
in predicting the maximum potential yield. 
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(3) Modeling of reflectance (back scattering) in the microwave region# 
This could be useful in the identification of small grains crops and 
the determination of surface soil moisture, plant moisture and the leaf 
area index. 

(4) Modeling of atmosphere sepreately as well as in tandem with the 
crop canopy reflectance. This obviously is very relevant since the 
eventual goal of the modeling effort is to use satellite based reflec- 
tance data to assess vegetation. 
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